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ABSTRACT: 

Chebyshev 's method has not got much attention in recent years compared to Newton 's 
method. Some of the papers which appeared recently discuss modification of Chebyshev 's method 
free from second order derivative. In this paper, we have modified Chebyshev's method by suitably 
approximating second order derivative using Taylor's Series. The proposed method requires 
evaluation of only three functions and still maintains cubic convergence. For a particular choice of a 
parameter in the method, we get fourth order convergence. Examples are provided to show the 
efficiency of the method with classical Chebyshev 's method and few other cubic convergent methods. 

KEYWORDS: Chebyshev's Method, Iterative Method, Non-linear equation, Second derivative free 
method, Cubic convergence. 



I. INTRODUCTION 

One of the most important and challenging problem in applied mathematics and engineering is to find 
an approximate solution of the nonlinear equation 

fix) = (1) 

where f:D cK-il for an open interval D is a scalar function. Newton's method is one of the famous 
iterative methods to solve (1). It is well known that it has quadratic convergence. The iterative formula of 
Newton's method is given by 

3fe=*»-£8.*=<U.2._ (2) 

By improving Newton's method, Chebyshev [8] derived a cubic convergent iterative method which 
requires computing of second derivative of the function (1). Without evaluating second derivative of f(x), cubic 
convergence was first established by the Arithmetic mean Newton's method [1, 16]. In [1, 4, 11, 14] third 
order accuracy is proved by using midpoint Newton's method. Modifications in the Newton's method using 
harmonic mean were suggested in [1, 7, 11, 14]. Nedzhibov [13] gave several classes of iterative methods 
using different quadrature rules for solving nonlinear equations. Different variants based on Simpson formula 
on Newton's theorem were proposed in [1, 2, 5, 18]. Zhou [17] gave a class of Newton's methods based on 
power means on the trapezoid formula. Tibor et al [12] gave third order accurate Newton's modification by 
using geometric mean. Jisheng et al [10] proposed a uniparametric Chebyshev type method free from second 
derivative. Herceg et al [6] presented a method for constructing new third order methods for solving (1) in 
which Halley's and Chebyshev's methods are special cases. Esmaeili et al [3] gave a uniparametric 
modification of Chebyshev's method. Different variants of Chebyshev's method with optimal order of 
convergence were proved in [15]. Jayakumar et al [9] modified Newton's method using the harmonic mean 
instead of arithmetic mean on the Simpson's formula. The methods presented in all the above papers require 
only the first derivative of f(x) and establishes third order convergence. 

In this paper, we propose a modification in the Chebyshev's method [8] similar to the one given in [10] 
which is free from second order derivative of f(x) by considering a suitable approximation of f "(x) using 
Taylor's Series. The proposed new method has the advantage of evaluating only the first derivative of f(x), less 
number of iterations and third order accuracy. With an additional condition on the parameter k, fourth order 
accuracy is obtained. In Section 2, we present some definitions related to our study. In Section 3, some known 
third order variants of Newton's method are discussed. Section 4 presents the new method and its analysis of 
convergence. Finally, Section 5 gives numerical results and discussion. 
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II. DEFINITIONS 

Definition 2.1 [16]: Let s B R, x m E fi, n = 0,1,2, . . . Then the sequence {.r n } is said to converge to a if 

li'ml*- — ctl = 0. If, in addition, there exist a constant c > 0, an integer n D > and p > such that for all 

n-si 

B> ) l-r^+i-al < ck rl -or| p , 

then {^ffi}is said to converge to a with order at least p. If p = 2 or 3, the convergence is said to be quadratic or 
cubic. 

Definition 2.2 [8]: The efficiency index of an iterative method in the sense of Traub is defined by the 
equation E' = p '" , where p is the order of the method and n is the total number of function and derivative 
evaluations at each step of the iteration. 

Definition 2.3 [16]: Let a be a root of the function (1) and suppose that * n _i., x v an d *n+L are three 
successive iterations closer to the root a. Then, the computational order of convergence (COC) p can be 

. , . . u c , _. iDli;j,;-i-iri:'i;j r .-i^i| 
approximated using the formula Q * —. ; r . 

III. KNOWN THIRD ORDER METHODS 

Let a be a simple root of a sufficiently differentiable function f(x~). Consider the numerical solution 
of /Of) = 0. From Newton's theorem, we have 
/Of) = f(x n :\ +£ f{t)dt (3) 

Arithmetic Mean Newton's method (AN): By using the trapezoidal rule in (3), we obtain 

£^^[/w + iifiivfc -^) +/oo] ( 4) 

From equations (3) and (4), a new approximation x a +i for x is obtained and for m = 1, we get the AN [1, 16] 

*- = ^-^Sw n = at2 - (5) 

Since equation (5) is implicit, we can overcome this implicit nature by the use of Newton's iterative step (2). 

Harmonic Mean Newton's method (HN): If we use the harmonic mean instead of arithmetic mean in (5), we 
obtain harmonic mean Newton's method [1,7, 11, 14]: 

*-i = *» - tfbfcVW ' " = °' lj 2 '"■ (6) 

where y„ is calculated from (2). 

Midpoint Newton's method (MN): If the integral in (3) is approximated using the midpoint integration rule 
instead of trapezoidal rule, we get the midpoint Newton's method [1,4, 11, 14]: 



-, n = 0,1,2,..., (7) 



/*(- — -) 

2 

where 34 is calculated from (2). 

IV. MODIFICATION OF CHEBYSHEV'S METHOD AND ITS ANALYSIS OF 

CONVERGENCE 

The classical Chebyshev's method (CCM) is given by 

Use Taylor's series to approximate/ OfjJ in (8) with — — , where >"„ is given by 

.Via - *n 
y = x -k — , k * 0, n = 0,1,2,... (9) 

^ n n 15 5 5 5 \ y 
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Thus equation (8) is converted into 

We can modify equation (10) by using (9) as below 



(10) 



»*-|l-i 



/(* n > 



/ 






/"Ok.) -/'<*„) 



fC«n3 






where the parameter k ■*■ . Finally, we simplify the above equation to obtain the new method 

{ i r /■(>.)- /'u„)~n /u.) 

x + I -1 + 1 = — I — , n = 0,1,2,... . 

2 * L /'»») J J /'(*») 



(11) 



Theorem 4.1 Assume that the function / :J3 <= K -» EL has a simple root aei?., where D is an open interval. If 

fix) is sufficiently differentiable in the interval D, then the methods defined by (1 1) converge cubically to k in 

a neighbourhood of a. for all values of the parameter k * . Moreover, if / "(# ) = and k = 2/3 then the 

method has fourth order convergence. 

Proof. Let a be a simple zero of function /Of] = 0, (That is, /(h) = and / '(«) * 0). 

Let Bf, = x n — a and d n = v n — x^ .. where 3$, is defined by (9). By expanding /(cO by Taylor's series 

about -Tf,. we obtain 



/(«) = /C*J - **f(xj+lf"(x»>4 - \fbJ4 + ofe*>. 

Taking into account /(a) = 0, we have from (12) 



(12) 



(13) 



Let us denote c 



;t*j 






j = 2,3,.... Dividing equation (13) by / Of,,), we get 



From (9) and (14), we have 

&% - y : v. - K * - ~^77 n \ - _feff n + kC-4 - kCi4 +" ^i4. + ®(4~}- 

Using Taylor's expansion for / (}■„] about *,,, we get 



(14) 



(15) 



(16) 



From (15) and (16), we get 

/' <yJ - /W = fbJ k [-e n + Cz4 ~ C 3 s| + C*4] + 

-2«7 4 + C 2 CJ4] + ^^ k*[-4 +3(.C 2 4 - C\4 + C^>] + ^^fc>* - 4C 2 fl|] + °( eJJ). 



^fc a [^ - 2C 2 4 + (C| + 2C,)4 



(17) 



Dividing (17) by / \x) and the result is multiplied by — to get 

Ik 
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ik (^^^f^) = - e "- C = + e - ( C ^ + l ftC - s - C ^^ (l + 3fe:) + ^fc = ^] + ^S [3fc (i^C 3 + 2C#) + 



C 2 C 4 (1 + 6k 2 ) + | fc J C 5 ]+ G(b-J) 



(18) 



From (14) and (18), we have 

Again from (14) and the above equation, we get 

-^S + 5[( £a ^)]^ = -" + 'j( 2 ^ + (T- 1 ^)-»S(v.P + l") + ^^.^'- 



Using the relation s n = x„ — a in (1 1), we have 

Finally, using (19) on (20), we obtain 

a^i = 4 [2Ci + C 3 (^ - l)] - 4 [ci + C 2 C 2 (2 +^) + C 4 (2fc J - 1>] + QOg) 



(19) 



(20) 



(21) 



From equation (21), we conclude that the methods defined by (1 1) are cubically convergent for all the values of 
k ^ 0. Also, taking / "(x ) = (i.e. C 2 = ) and k = 2/3 in equation (21), we get fourth order convergence. 

V. NUMERICAL EXAMPLES AND DISCUSSION 

In this section, we present the results of numerical calculations to compare the efficiency of the present 
method (MCM). We compare MCM with Newton's method (NM), classical Chebyshev method (CCM) and 
other cubic convergent Newton-type methods given in section III. Number of iterations and Computational 
Order of Convergence (COC) for different functions is given in Table I. Numerical computations are worked 
out in the MATLAB software with double precision accuracy and the results are presented below. 
Stopping criterion for the iterations has been taken as U n+ i_ — x n l + lff%+i)l < e, where s = 10 -14 . The 
following table gives functions that are taken as examples: 



ft*) 


Root (a) 


ft*) 


Root (a) 


fit*) 


= arctan(x) 





fit*) = in (* J + x + 1} 





/zM 


= e~*sinx + log(l + x 2 } 





f 5 (s) = x 2 sin(x) - cos(x) 


6.3083089552381 


£t*> 


--x 3 -9x 2 +28x-30 


3 


f 6 (x} = (x + 2)e-x + x 


-1.6878939988 



Discussion: In this work, a single parameter family of Chebyshev-type method free from second order 
derivatives for solving non-linear equation is presented. Third-order convergence is established for the 
proposed method for all values of the parameter k ^ 0. Also, we get fourth order convergence when k = 2/3 and 
C2 = which is displayed for the examples f L (x), f 2 (x"j and/ 3 (at). The efficiency index for MCM is E = 
1.4422 and for the value of k = 2/3 and C2 = 0, E = 1.5874 which are very good compared to Newton's method 
which has the efficiency E = 1.4142. This shows that MCM is preferable to Newton's method, especially when 
the computational cost of the first derivative is not more than that of the function itself. These proposed 
methods are also superior as there is no evaluation of second derivative compared to its classical predecessor, 
Chebyshev' s method. From numerical examples, we have shown that these new methods have very good 
practical utility. 
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Table I 


m 


Xo 


Number of iterations 


COC 


NM 


AN 


MN 


CC 

M 


MCM 
k=l 


MCM 
k=l/2 


MCM 

k=2/3 


NM 


AN 


MN 


CCM 


MCM 
k=l 


MCM 
k=l/2 


MCM 
k=2/3 


fi(x) 


1.3 


8 


6 


5 


8 


6 


NC 


NC 


3.0 


2.9 


3.0 


3.0 


2.9 


NC 


NC 


1 


6 


5 


5 


6 


5 


6 


5 


2.9 


2.9 


3.0 


3.0 


2.9 


3.3 


4.9 


0.5 


5 


4 


4 


5 


4 


4 


4 


2.9 


2.9 


3.1 


3.0 


2.8 


3.5 


4.9 


-1 


6 


5 


5 


6 


5 


6 


5 


2.9 


2.9 


3 


3.0 


2.9 


3.3 


4.9 


m 


1.3 


5 


4 


5 


5 


5 


5 


4 


3.0 


2.2 


2.9 


2.56 


2.9 


2.8 


4.0 


1 


5 


4 


4 


5 


4 


4 


4 


3.0 


2.6 


2.8 


2.79 


2.5 


3.0 


3.6 


0.5 


4 


4 


4 


4 


4 


4 


4 


2.4 


2.6 


2.7 


2.78 


2.6 


2.7 


3.8 


-1 


6 


5 


5 


5 


5 


5 


5 


3.0 


3.1 


3.0 


3.07 


3.1 


2.8 


4.4 


Mx) 


2 


6 


5 


5 


5 


5 


5 


5 


2.9 


2.9 


2.7 


2.52 


3.1 


3.0 


4.8 


2.5 


7 


6 


5 


6 


6 


6 


7 


2.9 


2.9 


2.9 


2.97 


3.0 


2.8 


4.8 


3.5 


8 


6 


6 


7 


7 


7 


6 


3.0 


2.9 


2.9 


1.51 


2.6 


1.5 


4.8 


1 


9 


7 


7 


7 


7 


7 


7 


2.9 


2.9 


2.9 


2.78 


3.0 


2.9 


4.8 


Mx) 


1.3 


6 


4 


5 


5 


4 


5 


4 


2.0 


1.8 


3.1 


2.9 


1.9 


2.9 


2.2 


1 


6 


4 


4 


5 


5 


4 


4 


2.0 


2.0 


2.6 


3.0 


2.9 


3.5 


2.3 


0.5 


6 


4 


4 


4 


4 


4 


4 


2.0 


2.5 


3.1 


2.5 


2.4 


2.7 


2.4 


-1 


8 


7 


6 


7 


7 


7 


7 


2.0 


3.0 


3.0 


2.98 


3.0 


3.0 


2.9 


fix) 


4 


6 


4 


4 


5 


4 


5 


4 


1.9 


2.4 


2.5 


2.9 


2.2 


2.9 


2.3 


1 


5 


4 


4 


4 


4 


4 


4 


1.9 


2.9 


2.9 


2.9 


2.9 


2.9 


2.9 


0.5 


7 


5 


5 


7 


9 


8 


7 


2.0 


3.0 


3.0 


3.0 


3.1 


2.9 


2.8 


-4 


6 


4 


5 


4 


5 


5 


5 


1.9 


2.1 


2.9 


2.05 


2.9 


2.9 


2.9 


Mx) 


-3 


8 


6 


6 


8 


6 


6 


6 


1.9 


2.8 


2.9 


1.99 


2.9 


2.9 


2.9 


-2 


6 


4 


4 


6 


5 


5 


5 


1.9 


2.8 


2.9 


1.99 


2.9 


2.9 


2.9 


-1 


9 


7 


6 


13 


15 


13 


D 


2.0 


3.0 


3.0 


1.99 


2.9 


2.9 


D 



NC - Not Convergent, D - Divergence, x - Initial point and COC - Computational Order of Convergence 
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